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I present an overview of some of the recent results on jets and jet-like correlation measurements 
from the Relativistic Heavy-Ion Collider (RHIC) at Brookhaven National Laboratory. 

Jets are produced in the initial hard scatterings of an event and can therefore be exploited as 
probes of the hot and dense medium produced in heavy-ion collisions. Previous RHIC results 
indicate that this medium, the Quark Gluon Plasma (sQGP), is strongly coupled, with partonic 
degrees of freedom. High pr colored partons passing through the sQGP are therefore believed to 
suffer energy loss via induced gluon radiation and elastic collisions, before exiting the medium and 
fragmenting in vacuum. Jet reconstruction and high px correlation studies allow us to investigate 
how the partons interact with the medium and how the medium responds to the partons moving 
through it. By comparing measurements from pp and d-Au to those in Au-Au collisions at ^snn = 
200 GeV we aim to disentangle cold nuclear matter effects from those of the hot and dense sQGP. 

I. INTRODUCTION 

Over the past decade the RHIC experiments have produced significant evidence that a Quark Gluon Plasma 
(sQGP) is being produced in ultra-relativistic heavy-ion collisions. This sQGP is strongly coupled and has 
partonic degrees of freedom. Hard probes are now being used to study how partons interact with the medium, 
and how the sQGP responds to energy deposited by these highly energetic partons as they travel through it. 
The results from Au-Au collisions are compared to those from pp and rf-Au collisions where no QGP is believed 
to be created. 

Jet quenching, the loss of energy of hard scattered partons to the medium, was first observed at RHIC via the 
single particles nuclear modification factor, Raa [I]- The nuclear modification factor is defined as the ratio of 
Au-Au pt spectrum normalized to the number of binary collisions (Nbi n ) to that of the pp data measured at the 
same collision energy. The Raa of high px particles in central Au-Au indicated a suppression of up to a factor 
of 5 PP, while photons, colorless objects, reveal an Raa =1 0- These results suggest that the suppression is 
due to partons interacting with, and hence loosing energy to, the hot and dense colored medium. Raa measures 
however have a couple of limitations. First it is not possible to infer the initial partonic energy from the final 
state hadron. Second surface emission leads to an inherent insensitivity to the medium's density; no matter 
how dense the medium those partons near the surface will always escape and be detected [3J. 

To attempt to alleviate both these issues the collaborations are now performing full jet reconstruction. An 
added advantage is that by studying the fragmentation patterns we can hope determine not only how much 
energy is lost from the initial parton but also how the energy is redistributed. Both STAR and PHENIX use 
jet finding algorithms from the Fast Jet package [I], in addition PHEINX uses a Gaussian filter code, details of 
which can be found here [5] and STAR has used a traditional mid-point cone algorithm [5] . 

Neither RHIC experiment has hadronic calorimetry included in their detector designs, hence the full jet energy 
is not directly accessible but is assessed by combining charged particle momenta via tracking devices with neutral 
energy measurements from electro-magnetic calorimetry. In addition the energy from long-lived neutral hadrons, 
such as the neutron and JC£ is missed. This leads to a significant difference in the jet energy reconstructed and 
that of the initial, so-called particle level jet. This difference, as well as the jet energy resolution, has to be 
evaluated and corrected for before final results can be presented. For instance, the detector performances have 
been evaluated via simulations by STAR and show that the jet energy resolution in pp data varies from 10-25%, 
for 40-10 GeV/c jets [ZJ. 

II. JETS IN PP AND D-AU COLLISIONS 

For pp collisions the reconstructed raw jet spectra reconstructed with the Anti-K-r, Kt and SISCone algorithms 
were the same within 10%, confirming that they have similar behaviours in this low multiplicity data, Fig. ref- 
Fig:JetCompare [7]. The inclusive jet, Fig. [2] and di-jet cross-sections, Fig. |3| have been measured by STAR 
using the increased statistics of the 2006 data |S]- A midpoint cone algorithm |J5 with a cone radius of 0.7, 
a split-merge fraction 0.5 and a seed energy of 0.5 GeV was used. When hadronization and underlying event 
uncertainties are included both sets of data are well described by NLO theory [TO] ■ 
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FIG. 1: The ratio of reconstructed jets for various FastJet algorithms, Kr/Anti-KT and SISCone/Anti-KT as a function 
of reconstructed jet pr- 
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FIG. 2: Left: The 2006 measured inclusive jet cross section for pp collisions at -^/sjviv=200 GeV. Right: (Data- 
Theory) /Theory. 



The measured fragmentation functions of both STAR [TT] and PHEINX [T2j agree, within errors, with 
PYTHIA simulations |13j, see Fig. [I] The PHENIX data use the Gaussian filter with a width of 0.3, and 
are fully corrected. While the data from STAR are not yet corrected to the particle level and are therefore com- 
pared to PYTHIA 6.410 [13], tuned to the CDF 1.96 TeV data (Tune A), predictions passed through STAR's 
simulations and reconstruction algorithms. The agreement with PYTHIA simulations even for R=0.7 suggests 
that there are only minor NLO contributions beyond those mimicked in the PYTHIA parton-shower calculations 
at RHIC energies. 

We also investigate jet production in d-Au collisions, where cold nuclear mater effects are expected to be 
present but no QGP formed. For example, the presence of the Au nucleus may induce additional initial and 
final state radiation, or result in scatterings of fragmentation particles as they escape the nucleus. The measured 
d-Au mid-rapidity jet nuclear modification factor, Rep , where peripheral (i-Au events are used instead of pp 
data, is shown from PHENIX in Fig. [5] left panel, for three different centrality bins. Here the Anti-Ky algorithm 
with R=0.3 was used. A slight modification of the jet spectrum is observed in ti-Au collisions, with the central 
d-Au jet cross-section showing the greatest suppression. These results are consistent with the tt° results and are 
likely an indication of cold nuclear matter effects such as modifications of the nuclear PDFs and/or energy loss 
in the cold matter. Since these effects may result in more subtle modifications than that of the overall jet yields, 
another way to probe for these cold nuclear matter effects is via di-jet correlations. Re-scatterings in the nucleus 
my result in a broadening of the di-jet A(p distribution. The mean transverse momentum of the fragmentation 
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FIG. 3: Left: The 2006 measured inclusive di-jet cross section for pp collisions at -^/Sjviv=200 GeV. Right: (Data- 
Theory) /Theory. 
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FIG. 4: Left: PHENIX Run-5 pp charged particle (electrons rejected) fragmentation function as a function a z (=phadron 
/pjet- A vertical scaling of e(p^ ) = 10 fe ,fc = 0, ... ,3 is applied. The shaded boxes indicate systematic uncertainties, 
and error bars indicate statistical uncertainties. From [12] , Right: Jet FF as a function £ from STAR using the Kt, 
Anti-Kr and SISCone algorithms for R=0.7. The curves are the PYHTIA predictions. 



products with respect to the jet axis, {jr), and the mean transverse momentum kick given to di-jet pair, (fcr), 
are two variables used for these investigations. The (jr) was measured via di-hadron correlations and found 
to be constant at «0.55 GeV/ c for all pt triggers measured and for both pp and <i-Au events at y / s]v r /v=200 
GeV [15] , However, as shown in Figure [5] right panel the (fcy) is systematically higher for the d-Au data than 
for the pp data for all px jet and pt trigger ranges measured. This suggests that while cold nuclear matter 
effects are small they still result in a minor deflection/broadening of partonic trajectories, the fragmentation 
appears to be unaffected. 

The underlying event (UE) is an important element of hadronic collisions and is defined as those particles not 
produced in the initial hard scatterings. Hence it includes beam-beam remnants, particles from initial and final 
state re-scatterings and those resulting from soft or semi-hard multi-parton interactions, pile-up is not included 
in the UE definition and must be removed. In pp events at RHIC the UE is small and is often neglected, however 
in d-Au collisions it becomes sizable. CDF initiated such an analysis [TS]. First the jets are reconstructed, next 
each event is split into four sections defined by their azimuthal angle with respect to the leading jet axis (A0). 
The range within | A</>|<60° is the lead jet region and an away jet area is designated for | A0|>12O°. This leaves 
two transverse sectors of 6O°<|A0|<12O° and 12O°<A0<-6O°. One is called the TransMax region and is the 
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FIG. 5: Left : Jet Rcp for d-Au collisions for three different centrality selections. Right: The measured (kr) for d-Au 
and pp data at ^/smn—^OO GeV from di-hadron correlations and di-jet measurements from STAR. Vertical bars show 
the statistical and systematic uncertainties, horizontal bars indicate the bin widths. 
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FIG. 6: Mean pr (left) and mean number of charged particles per unit n and <f> (right) in the transverse regions for pp 
and d-Au collisions. All data are from -^/sjvjv=200 GeV. 



transverse sector containing the largest charged particle multiplicity. The second sector is termed the TransMin 
region. Two sets analyses can then performed, a "leading" jet study, where at least one jet is found in STAR's 
acceptance, and a "back-to-back" study which is a sub-set of the "leading" jet collection. This sub-set of events 
has two (and only two) found jets with p^" ay: > et /p£ a : ' et >0.7 and |A0j et |>15O°, this selection suppresses hard 
initial and final state radiation of the scattered parton. The TransMax region has an enhanced probability of 
containing contributions from these hard initial and final state radiation components. Thus, by comparing the 
TransMax and TransMin regions in the "leading" and "back-to-back" sets we can extract information about 
the various components in the UE. The properties of the UE in both pp and d-Au events are being studied, 
this is the first time such an analysis has been undertaken for d-Au collisions. Since this study is preformed at 
mid-rapidity it is likely that there is little to no contribution from the beam-beam remnants. Both the number 
of particles in and the momentum distribution of the underlying event appear to be largely independent of 
the leading jet's pt in both pp and d-Au collisions, Fig. [6] The mean transverse momentum is similar for pp 
and d-Au events in both the TransMax and TransMin regions as can be seen in Fig. [6] left panel. Meanwhile 
the average number of charged particles per unit rj and <fi increases by ~factor 5 from pp and d-Au collisions, 
right panel of Fig. [6j This increase in particle production is only slightly less than N part scaling of the pp 
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data would predict, also shown in figure [6| All the results are the same within errors for the "leading" and 
"Back-to-Back" data sets [13 [18], which suggests that the hard scattered partons emit very small amounts of 
large angle initial/final state radiation at RHIC energies. 



III. JETS IN HEAVY-ION COLLISIONS 



The presence of jets in heavy-ion events is clearly evident in Fig.[7| despite the significant underlying event. 
The left panel shows a di-jet event in the PHENIX detector from a ~ 20% central Cu-Cu collision as found by 
the Gaussian filter algorithm [15]. The right panel shows a central Au-Au event in the STAR detector. Each 
grid cell indicates the summed px of the charged tracks reconstructed in the TPC and neutral energy recorded 
in the electro-magnetic calorimeter. 




FIG. 7: Di-jets in heavy-ion events at ^/s^ =200 GeV Left: PHENIX Run-5 Cu-Cu ~ 20% centrality. Charged tracks 
and photons are shown at the bottom by a lego plot. The distribution of the Gaussian filter output values of the event 
is shown at the top as a contour plot. The maxima in the filter density are reconstructed as jet axes, shown as red lines 
at the positions on the contour and Lego plots. Right: STAR central Au-Au event, the summed pt of charged tracks 
and electro-magnetic calorimeter towers per grid cell are shown. Clear di-jet peaks emerge above the background. 



In order to extract information regarding jets and the interactions of hard scattered partons with the sQGP 
it is essential to first understand the enormous background the jet is immersed in and its fluctuations. This 
background is predominantly formed from the soft underlying bulk particle production and is thus strongly 
dependent on the centrality of the collision. Schematically, assuming all the initial partonic energy is recovered 
by the jet finders the measured jet spectrum in Au-Au collisions is: 

do-AA do~AA _ „, , , M « 

— = — ®F(A,p T ) (1) 

apr dpr 

where F(A, px ) accounts for the background and its fluctuations, which are a function of the jet area and the 
Pt of the reconstructed jet. Initially it was assumed that these fluctuations could be accounted for by a simple 
Gaussian ansatz, however more detailed studies have shown this modeling to be insufficient. If the background 
is due to independently emitted particles then in a fixed area the number fluctuations are well described by a 
Poisson distribution, while those of the mean px result in a Gamma distribution, assuming a fixed number of 
particles, M [5T]- Therefore 



F(A,pr) = Poisson(M(A)) ® T(M(A)(p T )) (2) 

Such a modeling of F(A, px ) gives a good description of the the summed px in the random cones of area A 
on a toy simulation where particles with dN/d?7 = 650 are thrown with a T=290 MeV. The mean of the F(A, 
px) distribution is given by pA. p is the median{py e *' reco / J 4i} in the event and p 3 ^ t,reco is the reconstructed jet 
Px-- When however the FastJet jet finders are used, the description is not exact, suggesting that the jet finders 
clustering does not occur in a truly random fashion, as should be expected. To investigate the resilience of the 
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jet finding in heavy-ion events to these fluctuations we are using probe embedding into real Au-Au events. A 
particle, or jet, of a known transverse momentum, p^ nbed , is embedded into an event, the Anti-K-r algorithm 
run and the reconstructed jet containing the embedded probe identified. Then 

£~ ~J et ; reco A embed fo\ 

dp T = - p- A-p T (3) 

is calculated. The 5px distributions are then calculated over many events for different embedded objects. Single 
pions, PYTHIA jets, and qPYTHIA jets (where the fragmentation pattern is altered from that of vacuum 
fragmentation) [55] have been used as well as various p|Z™ bed . Figure pi left panel shows the resulting 5px 
distribution for a single 30 GeV/c pion embedded into a 0-20% central Au-Au event [23]. The right panel of 
Fig. [H] shows the distribution of the event-wise difference 

A5p T = 6 P Z-S$ t , (4) 

between S P t for a PYTHIA-generated jet with pt > 30 GeV/c probe and that of a pion with the same pt,Vi 
and (/). Similar results are seen when qPYHTIA jets are used and for lower px probes. This study reveals 
the Anti-KT response is insensitive to the fragmentation pattern of the probe, greater than 70% of the time 
| ASpx\ < 200McV/c. This is a crucial property for the jet finder used in heavy-ion jet quenching studies since 
we do not yet have a complete description of the fragmentation functions of partons which traverse the sQGP. 
Also importantly the Anti-Ky algorithm has been shown to respond in a predominantly geometric fashion and 
hence is fairly impervious to back-reaction effects 24 . 




FIG. 8: Left: Spr distribution for a single pion with px = 30 GeV/c embedded into a 0-20% central Au-Au event. Right: 
The event-by-event difference in 8pr for a PYTHIA embedded jet, px >30 GeV/c, and a single pion with same px, n 
and 4> as the jet. 

For unbiased jet reconstruction one would expect the jet Raa to be close to unity, possible deviations might 
occur due to initial state effects in the Au-Au collisions. However, the left panel of Fig. [9] shows that even for 
R=0.4 the jet Raa is likely below unity, due to the large systematic errors the results are just compatible with 
unity |25j . The jet Raa is however significantly above that of single hadrons withp^ >20 GeV/c (i?^™ 11 ~ 0.2). 
One can also observe that there are significant differences between the results for the Kx and Anti-K-r algorithms 
which is expected given their different responses to the heavy-ion background. The ratio of the number of 
reconstructed jets for R=0.2 compared to R=0.4 is less for Au-Au data than for pp right panel of Fig. [9] |25j . 
Taking these two results together we conclude that the jet algorithms do not recover as much of the original 
partonic energy in Au-Au events as the same algorithms and settings run on pp data. Further, Fig. [9] indicates 
that this is likely due to the fact that particles are emitted at larger cone angles in Au-Au events compared to 
pp events with the same jet energy, with considerable energy, even at higher jet pt outside of R=0.4. 

Di-jet coincidence rate measurements provide us with evidence that there is a path-length dependence to 
the partonic energy loss. In this analysis "trigger" jets are identified which have a reconstructed jet px > 20 
GeV/c when only particles with px > 2 GeV/c are considered by the Anti-K^ algorithm. These trigger jets 
also contain a barrel electro-magnetic calorimeter tower with E tow > 5.4 GeV/c. This high z fragmentation 
requirement biases the trigger jet to being preferentially emitted from the surface of the medium and/or to 
have only minimally interacted with the medium. Such a surface bias in turn maximizes the average distance 
traversed by the recoil jet through the medium. If partonic energy loss is dependent on the path-length through 
the medium, the recoil jet will therefore reveal a greater suppression than that observed for the unbiased jet 
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FIG. 9: Left: Jet Raa- Right: The ratio of the number jets found with R=0.2 to those found with R=0.4 for pp and 
Au-Au events as a function of jet px- 



population. The relative probability of reconstructing a di-jet pair in Au-Au is compared to that in pp is shown 
in Fig. 10 This relative probability is suppressed by an approximate factor of 5 26J, i.e. a much stronger rate 



of suppression than observed for the inclusive jets. This results supports the notion of a path length dependent 
energy loss term. 
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FIG. 10: Left: Di-jet reconstruction probability in Au-Au over pp for R=0.4 using the Anti-Ky algorithm. The curves 
indicate the systematic uncertainties in the estimation of the background fluctuations (solid) and of the pt of the trigger 
jet assuming a jet resolution of 25% (dotted). Right: Jet-hadron correlations for 0.5 < p T ssoc < 1 GeV/c for pp and 
Au-Au collisions. 



To investigate further the jet broadening and softening indicated by the studies mentioned above we turn 
to jet-hadron correlations. In this analysis a "trigger" jet (defined as in the di-jet analysis above) is used 
to determine the jet axis and the A<fi correlation of all charged particles in the event relative to this axis is 
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examined, for more details on this analysis see [28 . An example of such a correlation is shown in Fig 
right panel, again a softening and broadening of the distributions of particles from jets is indicated for low px 
associated. The per trigger A0 distributions for pp and Au-Au event, plotted as a function of the associated 
charged particle px , are summarized in Figures [TT] and |12| The Gaussian widths of the away-side correlations 
in pp and Au-Au are shown in Fig, |11| The Au-Au distributions are broader than those in pp for low pt 
associated particles, accompanied by an significant increase in the low px associated yields |28j . For high px 
associated particles the Au-Au recoil jet correlation width is equivalent to that of pp but there is a significant 
reduction in the particle yield. Re-scattering of the initial parton could also potentially cause such a broadening 
rather than a modification of the fragmentation. Therefore the di-jet A<f> distributions in pp and Au-Au data, 
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PYTHIA events, and PYTHIA jets embedded into Au-Au events were studied. The results are shown in right 



plot of figure 11 The distribution is broader for the Au-Au data, however much of this broadening can be 
attributed to de-resolution of the jet axis due to the large underlying event, as a similar broadening is also 
observed in the PYTHIA+Au-Au event data. A similar result has been reported by PHENIX, who show that 
the A(f> distributions of di-jet events in Cu-Cu collisions do not vary as a function of centrality [57] ■ The red 
curve in the left plot of Fig. 11 indicates the expected width of the away-side A(f> distribution if the Au-Au 
fragmentation was pplike, but with the jet axis direction smeared to reproduced the width of the A</> Au-Au 
di-jet data. Clearly such a smearing cannot fully explain the observed broadening, and it also does not explain 
the enhanced low px yields. 
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FIG. 11: Left: The Gaussian widths of the away-side correlations as a function of py ssoc . Right: Di-jet A<f) distributions 
for various data sets at % /Sjvjv=200 GeV collisions. Solid curves are Gaussian fits to the data. 



The integrated yield difference, (Daa — Yield aa x (p T ssoc ) — Yield pp x (p T ssoc )) of the near- and away-side 



correlations as a function of p T ssoc are plotted in Fig. 12 As expected, the "surface" bias of the trigger causes 



the near-side Daa to be consistent with zero for all p^ ssoc . This means that there is an approximate energy 
balance, and a similarity of the associated px particle distributions for Au-Au and pp data for the trigger jet. 
The away-side data, Fig. [12] right panel, reveals that the low px hadron enhancement in the Au-Au data is 
approximately matched by the high px associated particle suppression. This suggests that the broadening and 
softening observed in the away-side correlation data is indeed due to a modification of the partonic fragmentation 
and not from residual soft background particles. 
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FIG. 12: Jet-hadron Daa distributions for the near-side (left) and away-side (right). For details on the systematic 
uncertainty bands see [28| . 



To remove the surface bias of the trigger object introduced in the di-jet and jet-hadron analyses discussed 
above PHENIX have been investigating 7-hadron correlations. Direct photon-hadron correlations are an ideal 
channel for studying energy loss since direct photons do not interact via the strong force and hence traverse 
the sQGP unmodified. At leading order pQCD, direct photons are produced from a Compton scattering of 
q + g — > q + 7 or quark annihilation q + q — > g + 7. To conserve energy and momentum a matching recoil 
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jet is also produced. The energy of the photon can then be used as a proxy for the jets initial energy and the 
fragmentation function of the recoil jet can be calculated. 7-hadron A(f> correlations have been measured in 
both pp and Au-Au collisions |29) . an isolation cut is used around the trigger photon to reduce contamination 
from 7T° and fragmentation photons. A fragmentation function is then deduced for the recoil jet correlation at 
\A(j> — 7T I < 7r/2. The resulting distributions as a function of £ = — Iu(xe) where xe = p T adron cos{A(j)) /p^ 10 * 1 " 1 
are shown in Fig. [13] for pp and central Au-Au collisions. The preliminary Au-Au and published pp data are 
plotted and compared to the TASSO measurement [30] and a Modified Leading Logarithmic Approximation 
(MLLA) in medium prediction [31] . The TASSO data and MLLA curve have been arbitrarily scaled down by a 
factor of ten to account for the limited PHENIX 77 acceptance as in [32] The shape of the isolated pp data are 
in good agreement with the TASSO measurement of the quark fragmentation function from e + + e~ collisions. 
The Au-Au results show depletion at low £ and possible enhancement at high £ compared to the pp data. These 
results again indicate that the energy lost by high px partons reappears as soft hadrons correlated with the 
initial parton's path. 




FIG. 13: £ distributions from PHENIX 7-hadron correlations in Au-Au data (black circles) and pp data (open blue 
circles) compared to TASSO data (green triangles) and MLLA in medium prediction (red curve). 



IV. SUMMARY 



In conclusion, both PHENIX and STAR are making quantitative steps in our understanding of jet production 
and fragmentation in pp rf-Au, and Au-Au collisions at RHIC at y'sjvw = 200 GeV over a wide kinematic range. 

Both the jet and di-jet cross-section in pp collisions are well described by next-to-leading order calculations 
once hadronization and the effects of the underlying event are taken into account. PYTHIA simulations repro- 
duce the measured distributions of the fragmentation products even at large jet resolution parameters indications 
that NLO corrections, beyond those implemented in PYHTIA, are small. 

Jet production in <i-Au collisions is slightly suppressed, particularly in central d-Au events when compared 
to binary scaled pp or peripheral <i-Au data. Together with the kx and ]t measurements as a function of jet 
Pt this indicates that there are small cold nuclear matter effects present but that these do not affect the shape 
and distribution of the fragmentation particles produced. 

Underlying event measurements of the pp and d-Au data show no significant changes as a function of jet p^The 
mean number of charged particles produced in the transverse region approximately scales with the number of 
participants in the event, while the mean px of these particles remains constant between pp and d-Au data. 

Our understanding of the background, and most importantly its fluctuations, in heavy-ion events has signifi- 
cantly improved. The Gaussian ansatz of the fluctuations has been shown to be incorrect, they arc more closely 
reproduced from a folding of a Gamma function with a Poisson that depend on the jet area, multiplicity of the 
background and its mean px- It has also been shown that the Anti-K-r algorithm's response to the background 
and its fluctuations is largely independent of the fragmentation pattern of the jet. 

Using a jet resolution parameter of R=0.4 the measured jet cross-section in central Au-Au collisions does not 
binary scale compared to pp data, the jet Raa < 1. This reveals that the lost partonic energy is spread to radii 
beyond R=0.4. Further the Au-Au R=0.2/R=0.4 ratio as a function of jet pr is lower than that measured in 
pp showing this broadening is there for all radii. The di-jet reconstruction probability in Au-Au collisions is 
suppressed as would be expected if the partonic energy loss is pathlength dependent. 
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Both jet-hadron and direct photon-hadron correlations indicate an enhanced production of hadrons a low px 
compared to pp baseline measurements which appear to compensate the suppressed particle production at high 
Pt- However, the high p^ associated particles while suppressed in number compared to pp data do not reveal 
any broadening in the jet-hadron correlations. This is in agreement with a scenario where the scatted parton 
loses energy in the sQGP but then fragments outside of the medium as it would in vacuum albeit with a reduced 
energy. The di-jet A(f> distributions indicate no obvious deflection of the parton's path although significant path 
length dependent energy is lost as it traverses the medium. 

All of the measurements made in heavy-ion collisions show that energy lost by high Q 2 scattered partons 
re-appears as soft particle production, with properties similar to that of the bulk, that remains largely correlated 
to the jet axis. 
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